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INERTIAL  MECHANISM  OF  SETTLING  OF  COARSELY  DISPERSED 
AEROSOL  ON  TERRESTRIAL  VEGETATION 


The  modern  theory  of  diffusion  of  a  heavy  additive  considers  only 
gravitational  settling,  i.e.,  it  actually  deals  with  the  settling  of 
particles  on  a  smooth  horizontal  surface.  The  fact  is  that  there  is 
really  no  ideally  flat  pl-ace  on  the  Earth's  surface;  the  latter  is  always 
rough  and  a  large  portion  of  it  is  covered  by  vegetation. 

The  settling  of  fine  particles  on  a  rough  surface  in  the  presence  of 
a  wind  cannot  be  governed  by  the  force  of  gravity  alone.  The  rate  w  of 
gravitational  settling  of  the  particles  is  usually  tens  and  hundreds  of 
times  less  than  the  wind  speed  u;  therefore,  the  particles  do  not  fall  to 
Earth  vertically  at  all,  but  move  above  the  ground  on  what  is  usually  a 
sloping  trajectory.  Along  with  the  gravitational  settling,  there  must  also 
be  a  settling  of  the  particles  on  plants  and  other  obstacles  under  the 
influence  of  the  forces  of  inertia;  under  certain  conditions,  it  can  also 
result  from  the  action  of  electrostatic  forces,  radiometric  forces  of 
Brownian  motion,  etc.  In  the  case  of  coarsely  dispersed  aerosols,  we  are 
interested  primarily  in  the  inertial  mechanism  of  settling,  in  addition 
to  the  gravitational  mechanism. 

This  can  be  supported  by  experimental  data.  In  an  experimental  study 
of  the  propagation  of  liquid,  coarsely  dispersed  aerosols  [1],  20  wind- 
vane  collectors  were  set  out  on  the  experimental  area,  at  different  dis¬ 
tances  from  the  aerosol  source,  to  evaluate  the  role  of  inertial  settling. 
Each  of  these  collectors  was  equipped  with  a  horizontal  glass  plate  and  a 
vertical  glass  slide.  The  wind  vane  turned  constantly  under  the  influence 
of  the  wind  and  the  vertical  slide  always  faced  into  the  wind.  A  horizontal 
glass  slide  was  placed  on  the  ground  beside  each  collector.  It  was  assumed 
that  above  the  horizontal  slide  in  the  wind  vane,  as  well  as  above  the  glass 
slide  located  on  the  ground,  a  laminar  boundary  layer  is  formed,  from  which 
the  drops  fall  onto  the  surface  of  the  glass  by  gravitation  with  a  velocity 
w,  and  that  the  deposition  on  the  windward  side  of  the  vertical  slide  occurs 
under  the  influence  of  inertia.  The  slides  in  the  wind-vane  collector  were 
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located  near  the  upper  limit  of  a  relatively  sparse  growth  of  vegetation, 
with  an  average  plant  height  h  -  30  cm.  for  each  fraction,  we  de  :ermined 
the  average  value  of  the  ratio  nvf for  20  points,  on  the  basis  of  the 

results  of  a  microscopic  analysis  of  the  slides;  in  this  ratio  is  the 

average  number  of  drops  of  a  given  fraction  which  land  on  a  unit  area  of 
the  windward  side  of  the  vertical  slide,  while  is  the  same  value  for 

the  top  of  the  horizontal  slide. 

The  results  are  shown  in  Table  1  in  the  form  of  the  ratio  =  fid) , 

where  d  is  the  average  diameter  of  the  drops  of  a  given  fraction. 


Table  1 

Average  Values  of  «y/n^  for  Various  Aerosol  Fraction,  Obtained  in  Experiments 
[1]  With  the  Aid  of  Wind  Vane  Collectors 


Number 
of  Ex- 

Di ameter 

d  of 

drops  ,  i n  mi  crons 

^  o 
a:  <u 

periment 

m 

13 

» 

1 

42 

5$ 

1 

!  75 

92 

103 

133 

167 

217 

s 

— '  10 

2 

o.si ! 

2.S 

4,3 

4,4 

4,5 

4,3 

4.1 

3.6 

2,5 

2,1 

2.2 

1.8 

3 

4.6 

3.3 

7,1 

7,5 

5.0 

3.2 

2,4 

3,0 

3.1  | 

2.9 

4.1 

3,6 

7 

Oeposi tion 
rate  of 

2,6 

5,2 

6.1 

5,6 

5,5 

4,4 

3.2 

2.4 

3,5 

1.2 

0,76 

1.6 

drops,  w, 
m/sec 

0,00550 

0,020$ 
1  ! 

0,0320 

0,112 

j  0,173 

0,234 

0,33$ 

0,467 

0,620 

0,953 

0,25 

— 

Note:  u[h)  is  the  mean  wind  speed  at  the  upper  limit  of  the  vegetation  with 
h  =  30  cm. 

It  is  clear  from  Table  1  that  the  values  of  n  /n,  >>  1,  i.e.,  on  the 

v  n 

vertical  slides  (which  can  be  viewed  as  crude  models  of  leaves),  the  density 
of  the  dr^ps  of  a  given  type  was  several  times  greater  than  on  the  horizontal 
ones . 


These  results  are  a  direct  indication  that  inertial  settling  of  aerosols 
on  vegetation  may  under  certain  conditions  prevaii  over  gravitational 
settling. 
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Having  adopted  this  point  of  view,  let  us  consider  the  propagation  of 
a  monodisperse  aerosol,  formed  by  a  constant  linear  source. 

Above  a  covering  of  vegetaion  with  height  h,  i.e.,  in  the  region 
0<x<‘x>,h<z<,x>{  with  source  coordinates  x  =  0,  z  -  H) ,  the  propagation 
of  the  aerosol  is  determined  by  the  equation  of  convective  diffusion: 


u(z) 


dq  (x,  z) 
dx 


dq  (x,  z) 
dz 


(1) 


where  q  is  the  concentration  of  the  additive  and  k  is  the  coefficient  of 
turbulent  diffusion. 

Turbulent  transfer  of  the  additive  takes  place  both  within  the  cover  of 
vegetation  and  above  it.  The  particles  settle  on  the  plants  by  inertia  and 
under  the  influence  of  gravity.  The  corresponding  decrease  of  additive  within 
the  cover  of  vegetation  can  be  thought  of  as  a  current  which  is  constantly 
flowing  through  the  plant  layer. 

The  inertial  current  of  particles  against  the  system  of  obstacles  around 
which  the  aerosol  flows  is 


£0=apu<7, 


(2) 


where  a  =  f(stk)  is  the  coefficient  of  entrapment. 
Stokes'  criterion  is 


slk  — 


(3) 


where  p 


P 


is  the  density  of  the  substance  of  the  particles,  d  is  their  diameter, 


u  is  the  viscosity  of  the  air,  x  is  the  characteristic  size  of  the  obstacles, 
a  1 

6  is  the  specific  (i.e.,  relative  to  a  unit  volume)  area  of  the  projection  of 
the  obstacles  on  a  plane,  normal  to  the  speed  vector  u. 


The  gravitational  current  of  the  additive  against  the  obstacles  is 


gff=afaq. 


(4) 
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where  B ^  is  the  specific  area  of  the  horizontal  projection  of  the  obstacles. 


In  accordance  with  this,  the  equation  of  turbulent  diffusion  for  the 
region  inside  the  cover  of  vegetation  (s  ^  h )  is  written  as  follows:1 


^  ^  -  -*-[*  « 

~  (a(~)?  (*)  u  (z)  -f  S h(z)  v  ]  q  ( x ,  z). 


(5) 


We  must  add  to  Equations  (1)  and  (5),  the  boundary  conditions  which 
determine  the  material  balance  of  the  additive,  the  conditions  on  the  ground 
(z  =  0) ,  both  at  infinity  and  contained  in  the  plane  z  -  h. 


To  solve  this  t..o-layer  problem,  we  must  know  the  microclimatic  and 
geometrical  characteristics  of  a  given  cover  of  vegetation:  a(s)  =  / 

[st/c(z0],  0(2)  =  f 2 (s) ,  8?i(sj  =  /3(s),  u(z)  =  f4(z,  Ri),  k(z)  =  f(z,  Ri)  in 

tn”*  region  0  <  z  <  h,  where  Ri  is  Richardson's  criterion.  At  the  present 
time,  these  characteristics  are  still  unknown2,  so  that  a  solution  of  the 
two-layer  problem  does  not  appear  to  be  possible  as  yet.  Other  approaches 
to  a  solution  must  be  sought. 


With  this  goal  in  mind,  let  us  consider  the  material  balance  of  the 
additive  for  an  element  of  the  cover  of  vegetation,  whose  length  is  dx,  whose 
width  is  equal  to  unity,  and  whose  height  is  h. 

With  a  stationary  regime  (continuous  source) ,  the  equation  of  material 
balance  of  the  additive  in  this  element  is  written  thus: 


J <?(-*.  z)atldS  z)dS  +  J  p(x,  z)dV «0, 

5  s  ■/ 


(6) 


where  z<^(a)  is  the  component  of  the  rate  of  turbulent  transfer  which  is  normal 

to  the  surface  Sj  H(XjZ)  is  the  component  of  the  diffusion  current  of  the 
additive  which  is  normal  to  the  surface  S,  and 


" p(x ,  z)  [o  (z]  i  ( 2 )  a  (z)  (z)  v]q  (x,z) 


1  Settling  under  the  influence  of  other  forces  can  be  considered  similarly. 

2  There  are  several  papers  dealing  with  the  wind  profile  within  the  cover  of 
vegetation;  in  particular,  data  for  corn  stalks  is  given  in  [2]. 
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is  the  rate  of  decrease  of  the  additive  within  the  cover  of  vegetation  due 
to  inertial  and  gravitational  settling;  S  and  V  are  the  surface  and  volume 
of  the  element,  respectively. 


Let  us  carry  out  the  integration,  averaging  over  z  and  disregarding  the 
small  higher  orders: 


|  <7  (x,  z)  u„  (c)  dS  =•=  jw  \q  (x,  h)  ~q(x,  0)1  —  huJl 
6  v  1  OX 

z)dS^k{h)JlixJildx 


oq  }  , 

w)dx> 


(considering  k( 0)  =0  when  disregarding  horizontal  diffusion): 

J  P  {x,  z)  dV  ~  —  J  (a  (z)  3  (c)  u  (z)  -f  3r  ( z )  w\  q  (x,  z)  dxdz  = 
==  —  (a  ^  tt  h  v  ^  jp/t)  q  dx. 

Substituting  these  expressions  into  (6),  we  will  have 

®  I<7  (x,  h)  -  q  (x,  0)]  -hu—L^-k  (A)  h)  - 
—  (a  g  uli  -f  Jffoh)  q  —  0. 


Let  us  replace  the  values  of  q,  3  /3  ,  and  u ,  averaged  over  s,  by  their 

VC 

values  at  the  upper  limit  of  the  cover  of  vegetation.  We  will  likewise  change 
the  values  a,  0,  and  g^ ,  calling  these  new  values  a,  g,  and  g^  .  Carrying  out 
the  indicated  substitutions  in  (7J,  we  will  have 


Oq  (x,  h ) 
Oz 


ha  ( h )  k\k2  Oq  (x,  h) 


q(x,  A)  + 


5 


where 


dq 

A  _  “  «. 

"'“ill*)1  K'2  ~~ T<fy  (a-TaT  • 

<?AT 

Disregarding  the  horizontal  gradients  q  in  relation  to  the  vertical  ones, 
we  finally  obtain  the  following  boundary  conditions  at  ti.e  upper  limit  of  the 
cover  of  vegetation: 


dq  lx,  A) 
dz 


aq  [x,  h), 


where 

_ (A)  +  ZJ  (h}h-  1+5) 

a= - - 

f_  q{x,  0) 

.  q  (x.  h)  • 


(8) 


(9) 


Hence,  the  problem  of  the  propagation  and  settling  of  an  additive  with 
consideration  of  both  gravitational  and  inertial  settling  leads  to  a  solution 
of  the  aquation  of  turbulent  diffusion  (1)  with  boundary  condition  (8)  on  the 
subjacent  surface3  for  the  value,  a  determined  by  Equation  (9). 

If  inertial  settling  prevails  over  gravitational  settling  (finfe  particles 
in  a  strong  wind,  tall  vegetation  through  which  the  wind  blows  readily),  then 
o8 u  (h)  »  B^w,  K  <  1,  and  we  can  use  the  following  expression  for  a: 


haiu  (A) 
A  (A) 


(10) 


Having  solved  the  diffusion  equation  with  a  new  boundary  condition,  we 
can  show  that  consideration  of  inertial  settling  by  the  root  method  changes 
the  picture  of  the  process. 

The  density  of  a  deposit  of  the  additive  on  the  ground  below  the  vegetation, 
produced  by  a  constant  source  during  1  second  (the  gravitational  flow  of  the 
additive  to  the  ground  below  the  plants) ,  is 

go=Zaq(x,h),  (11) 


3  Boundary  condition  (8)  was  proposed  in  its  general  form  by  Monin  [3] .  We 
gave  it  a  concrete  form  in  applying  it  to  the  inertial  settling  mechanism. 
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where,  according  to  (9) ,  5  -  "f, j'» 

The  density  of  a  deposit  of  the  additive  on  the  vegetation,  relative  to 
a  unit  area  of  ground,  produced  by  a  constant  source  during  1  second  (the  sum 
of  the  inertial  and  gravitational  flows  of  the  additive  onto  the  vegetation), 
is 


S,1ssSq  +gj~[kha+w (1  — 1)]<7  (x, k) . 


(12) 


For  practical  calculations,  it  is  necessary  to  know  how  to  determine  the 
new  dimensionless  characteristic  of  the  vegetation  which  enters  into  (9) , 
the  effective  coefficient  of  entrapment,  <y$h  =  f(stk). 

It  is  easy  to  see  that  the  direct  experimental  determination  of  this 
characteristic  for  an  actual  growth  of  vegetation,  for  example  by  direct 
measurement  of  the  additive  deposit  on  the  leaves,  involves  considerable 
difficulty.  Therefore,  let  us  use  the  indirect  method  of  determining  a$h 
on  the  basis  of  the  material  balance  of  the  additive: 

G~Go-r  Gh , 


where  G  is  the  output  of  the  source, 

a> 

G0  =  w?|^(x,  h)dx 

o 

is  the  gravitational  flow  of  the  additive  onto  the  ground  beneath  the  vege¬ 
tation,  and 


Gh  =  h  [<x?«  (//)  -fp?p]  f  q  (x,  h)  dx 
6 

is  the  sum  of  the  inertial  and  gravitational  flows  of  the  additive  onto  the 
vegetation.  Then 


Gh  __  G  —  Gq  _  a'iftu  (h)  ~  $,iyrh 

Go  Go  Zw  ~ 


1 


whence 


aVi  = 


“(*)'[*(  Go'  I)'~3//i] 


or  in  the  case  of  predominance  of  inertial  settling, 


(13) 


The  values  which  enter  into  the  righc-hand  side  of  (13)  can  be  determined 
by  means  of  field  tests  of  the  same  type  as  described  in  [1],  while  the  value 
w  is  calculated  in  the  usual  manner  [6] .  In  particular,  the  values  were 
determined  in  the  experiments  described  in  [1]  with  the  aid  of  wind-vane 


n. 


collectors  as  described  above.  It  was  assumed  that  £ 


where  n ,  is  the 
h 


average  number  of  drops  of  a  given  fraction  which  fall  on  a  unit  area  of  the 
horizontal  slide  in  the  collector;  n,  is  the  same  value  for  the  horizontal 

0 

slide  on  the  ground.  The  value  G  can  be  determined  in  field  tests  only  for  an 

L 

area  of  specified  dimensions,  i.e.,  it  is  assumed  that  tog  j  i]  (x,  h)dil^ 

oo  0 

»  j  <7 (*i A) dx,  where  L  is  the  width  of  the  area.  For  the  experiments 
o 


in  [1],  which  were  carried  out  with  the  source  at  a  height  of  1.6  meters  and 
an  area  500  meters  wide,  this  assumption  evidently  was  justified  for  all 
aei-sol  fractions  except  perhaps  the  very  finest  ones,  since  the  degree  of 
settling  e,  calculated  with  the  aid  of  the  Rounds  formula  (i.e.,  without 
considering  inertial  settling)  was  close  to  100%  for  the  majority  of  fractions; 
in  Experiment  No.  1,  for  example,  e  =  68.4%  for  fraction  d  =  50  microns, 

96.0%  for  d  =  83  microns,  etc. 

This  method  of  determining  the  characteristic  a&h  opens  up  the  possibility 
of  using  the  inertial  theory  for  practical  calculations. 
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Table  2 


Cover i ng  of 
Vegetation 

Rougnnes! 
Parametei 
z0  cm.  1 

5 

h  CM 

<x\h 

X  CM 

Sparse  grass,  3-10  cm 
high  [7]. 

0,52 

1,0 

5 

0,0015 

. 

Sparse  grass,  5-15  cm 
high  [7], 

0,66 

1,0 

8 

0,006 

_ 

Wi Id  plants ,  1959  [ 1 ] 

!  1,8 

1.0 

20 

0, 14 

0,20 

Wi  Id  pi  ants ,  I960. 

2,7 

1.0 

30 

0,055 

0,38 

Winter  wheat  in 

- - 

0,38 

60 

0,20 

0,30 

blossom,  50-80 
plants/m2 

Table  2  lists  the  values  for  the  parameters  5,  h  and  a &h  which  were 
found  in  field  tests  [1,  7]  and  in  additional  experiments  conducted  in  1%1 
on  wheat  sprouts.  As  we  know,  the  coefficient  of  entrapment  a  =  f(stk) 
approaches  unity  asymptotically  at  stk  ->  “.  Table  2  shows  values  of  a gh 
which  have  been  assumed  constant  for  the  sake  of  approximation  and  Correspond 
to  values  of  stk  >  8;  the  value  of  stk  was  calculated  by  (3) ,  while  the  values 
used  for  the  characteristic  size  of  the  vegetation  x  are  shown  in  Table  2. 

With  stk  <  8,  i.e.,  for  very  fine  drops  or  a  weak  wind,  the  values  of  cxgh 
must  be  made  smaller  than  those  shown  in  the  Table,  for  example,  by  multiplying 
the  values  in  the  Table  by  the  factor  stk/ 8. 

From  the  standpoint  of  practical  calculations,  an  important  advantage  of 
tne  new  method  (especially  for  agricultural  applications)  is  the  possibility 
of  calculating  separately  the  flow  of  the  additive  to  the  ground  below  the 
plants  (loss  of  chemicals)  and  the  flow  of  the  additive  directly  onto  the 
vegetation  (chemicals  usefully  employed) . 

For  practical 'approximate  calculations,  we  can  use  the  known  solutions 
of  the  diffusion  equation  (1)  q^Cx,  y,  0),  obtained  without  consideration  of 

inertial  settling  (for  example,  in  the  case  of  a  direct  linear  source,  the 
Rounds  formula  [4]  can  be  used),  but  the  additive  flow  must  be  divided  into 
the  two  parts  of  interest  to  us,  using  relationships  (1)  and  (2): 

Sh  + So  ^  ®<7o  (*.  0) h  [o$u  (A)  fyp  ]  q  (x,  h)  + 

+  ®  C*.  0)  =  { h  (ji)  -f  fyj®]  T  s® }  q  (*.  h), 


from  which  we  obtain  the  actual  concentration  of  the  additive  at  the  upper 
limit  of  the  covering  of  vegetation: 


q(x,  h) 


<7o  {x,  0) 


P 


ct3 u  (h) 


w 


ft  m9/m2 


Figure  1.  Density  values  for  deposits  of  an  additive  on 
the  ground  beneath  vegetation,  g q  =  f{x).  Experiment 

No.  1  [ 1 ] ,  d  =  25  mm. 

1  -  Experimental  data;  2  -  Exact  solution  of  diffusion 
equation  without  consideration  of  inertial  settling 
(Rounds  formula);  3  -  Rounds  formula  with  "inertial" 
coefficfient;  4  -  According  to  (15). 
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The  flow  of  additive  onto  the  vegetation  is 
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and  the  flow  of  additive  onto  the  ground  beneath  the  vegetation  is 
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(15) 


For  the  sake  of  an  example,  we  have  plotted  in  Figures  1  and  2  the 
results  of  determining  gQ  by  the  Rounds  formula  and  by  the  approximate 

formula  (15),  together  with  a  determination  of  <?Q(a:,  0)  by  the  Rounds 

formula.  The  following  values  for  the  parameters  were  used  in  the  calcu¬ 
lations:  (1)  Figure  1,  [1],  Experiment  No.  1,  fractions  17-33  microns, 
d  =  25  microns,  G  =  802  mg/m,  =  1.6  m,  7i  =  0.2  m,  k  =  0.166  m/secw 

w  =  0.0208  m/sec,  q  -  0.281,  u~  =  4.12  m/sec,  £  =  1.0,  a Bh  =  0.0486, 

y  2 

u{h )  =  2.62  m/sec,  k{h )  =  0.0332  m  /sec,  8^  =  0.  (2)  Figure  2,  [1], 

Experiment  No.  3,  fractions  67-82  microns,  d  =  75  microns,  G  =  3680  mg/m, 

V  -  1.6  m,  h  =  0.3  m,  k  -  0.234  m/sec,  w  =  0.173  m/sec,  q  =  0.312,  un  =  5.17  m/sec, 
£  =  1.0,  a &h  =  0.055,  u(h)  =  3.59  m/sec,  k(h)  =  0.0702  mr/sec,  8^  =  0. 

We  have  also  listed  here  the  corresponding  experimental  values  for  g q, 
obtained  in  experiments  [1]  with  liquid  aerosols. 

Obviously,  the  approximate  solution  is  significantly  closer  to  the 
experimental  data  tjian  is  Rounds'  solution. 

The  high  experimental  values  of  gQ  near  the  source  (x<  15  m)  can 

probably  be  explained  by  the  influence  of  the  aerodynamic  shadow  of  the 
automatic  device  on  which  the  aerosol  generator  was  mounted  [1]. 
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Figure  2.  Density  values  for  admixture  deposits  on  the 
ground  beneath  vegetation,  g ^  =  f{x) .  Experiment  No. 

3  [1],  d  =  75  microns.  For  key  to  symbols,  see  Figure 
1  . 


The  degree  of  settling  of  a  given  fraction  within  the  limits  of  the 
experimental  area,  determined  experimentally  (Figure  1),  was  e  =  4.18%; 
the  corresponding  values  according  to  approximate  inertial  formulas  are 
3.04%,  while  Rounds'  formula  gives  21.7%.  For  the  experiment  (Figure  2), 
the  experimental  value  of  e  =  36.6%,  while  the  calculated  values  are  38.9 
and  83,2%,  respectively,  i.e.,  when  inertial  settling  is  taken  into 
account,  the  calculated  values  of  c  are  much  closer  to  the  measured  values 
than  those  obtained  by  the  Rounds  formula. 

Moreover,  for  Figure  1  9^/Qq  =  6.73,  for  Figure  2  -  1.14,  i.e., 

as  in  the  experiments  with  wind-vane  collectors,  the  calculated  inertial 
flow  of  the  additive  onto  the  vegetation  exceeds  the  gravitational  flow 
for  small  (d  =  25  microns)  and  medium  (d  =  75  microns)  drops. 

Conclusions 

Direct  experiments  have  shown  that  inertial  settling  of  coarsely 
dispersed  aerosols  on  a  cover  of  vegetation  can  predominate  over  gravita¬ 
tional  settling  under  certain  conditions.  A  theory  is  proposed  for  the 
propagation  and  settling  of  aerosols,  which  takes  into  account  both  the 
gravitational  and  inertial  settling  of  particles  on  vegetation  covering  the 
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ground.  This  inertial  theory  makes  it  necessary  to  determine  a  new 
characteristic  of  such  cover,  the  effective  coefficient  of  entrapment. 

A  method  is  proposed  for  determining  this  characteristic  on  the  oasis 
of  field  test  results. 

A  comparison  with  experimental  data  shows  that  the  inertial  theory 
which  is  proposed  agrees  with  experimental  data  much  better  than  the 
generally  accepted  (gravitational)  theory;  in  addition,  the  inertial 
theory  makes  it  possible  to  determine  the  flow  of  the  additive  which  is 
deposited  directly  on  the  vegetation,  and  to  determine  separately  the 
flow  of  the  additive  onto  the  ground  beneath  the  vegetation,  which  is 
very  important  for  agricultural  purposes. 
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